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MODE COMPETITION IN THE QUASI-OPTICAL GYROTRON 


I. INTRODUCTION 

The quasi-optical gyrotron 1- ^ (QOG) is a promising millimeter and 

q 

submillimeter radiation source that has already demonstrated power 
levels < ISO kV vith pulse durations up to 13 ps at electronic 
efficiencies < 14%. Potential applications, ranging from heating of 
fusion plasmas to short wavelength radars, require stable CV operation at 
high power levels within a narrow frequency band. Because of the strong 
fields developed inside the cavity, the high density of longitudinal 

resonator modes, and the relatively long time of operation, nonlinear 

8 9 

interaction among many modes will inevitably occur. * In fact, it is 
the nonlinear coupling among cavity modes* 1 that will determine the 
existence of a final steady state, its accessibility from the initial 
small signal phase, and whether it is a single or a multimode state. 

In short, the evolution of the cavity fields can be described as 
follows. The injected electron beam initially excites all the cavity 
eigenmodes characterized by a start-up current below the beam current. 
At small amplitude each mode grows exponentially in time with its linear 
growth rate, unaffected by the presence of other modes. Soon one or more 
modes reach finite amplitude and start in**r acting through the induced 
modifications in the distribution of the .1 ctron beam. Vhen this 
happens a single mode may nonlinearly suppress all other unstable modes 
and eventually dominate. Equally well, it may destabilize modes that are 
linearly stable. In some cases the fastest linearly growing mode may be 
overtaken by another, linearly slower mode, causing mode switching. The 
final steady state, if one exists, may involve more than one large 
amplitude modes. The point to be made here is that the existence of a 
final steady state, be it single or multimode, is determined nonlinearly 
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and cannot be predicted from the linear behavior. 

Finding the regimes in parameter space associated with the desired 
type of final equilibrium proves too costly for direct numerical 
simulation, given that the performance of the QOG depends on four 
independent parameters. It should be stressed that the final equilibrium 
cannot alvays be obtained by examining the linear stability of every 
possible single mode at saturation.^ Even if a saturated single mode is 
found stable to small perturbations by other modes, there is no guarantee 
that the system will eventually evolve to this state. In general, more 
than one stable equilibrium exists for a given set of modes, however, 
only one of them is accessible from the appropriate initial conditions. 
In the case that the final steady state is dominated by a single mode, it 
is not necessarily the mode with the largest linear growth rate. 
Consequently, all participating modes must be treated on equal basis 
until (and if) one dominates. Mode selection and accessibility are 
inherently nonlinear processes. 

In this work an analytic model of the multimode dynamics is 
developed by expanding the nonlinear current in powers of the radiation 
amplitude. A set of equations is then obtained for the slow time scale 
evolution of the wave amplitudes and phases. The strength of the 
nonlinear interactions enters through coupling coefficients of known 
analytic dependence on the gyrotron parameters. Numerical simulations of 
particle trajectories are then replaced by a set of first order ordinary 
differential equations, one for each participating mode. Not only is the 
computation time reduced by orders of magnitude but, equally important, 
the final equilibria can be predicted analytically from the coupling 
coefficients. In this paper we will mainly discuss cases with up to 
three interacting modes, coupled up to fifth order in the wave amplitude. 
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The coupled equations obtained here can be applied to the description of 
12 

phase-locked operation of the QOG through the injection of a small 

external signal. They can also describe phase locking through the use of 

13 

a prebunching resonator, via certain modifications resulting from the 
use of a prebunched injected electron distribution in place of a uniform 
one. 
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II. GENERAL FORMALISM 


The configuration for the QOG is shown in Fig. 1. The cavity 
fields are expressed through the vector potential A(r,t), 


E = -cft’ B=7xA ’ (1) 

where A is a superposition of eigenmodes of the Fabry-Perot resonator, 

1 - -i» t 

A(r,t) - - £ A (t) U <r) e m e + cc , (2a) 

2 m y 

V r > - “p q ™ sk n>* ■ <2b> 

*m M m 

In the orthogonal coordinates (x,y,z) with x along the resonator axis and 
z along the electron beam axis, the transverse profiles u p qC r ) are given 
by the Hermite-Gaussian functions 


W 1 ' ‘ V (y) V 2) ’ 


h „ <*> ■ "p (iffy) ex p (—2~) ■ 

*m *m v ' w (x) 


(3) 


2 2 1/2 

v(x) - V (1+x /b m ) is the radiation spot size, V is the radiation 
2 

waist, b « k V /2 is the Rayleigh length and H (z) are the Hermite 
mm p 

polynomials. The wavenumber k m is equal to oWc, where the cavity 

frequencies « take the values 
m 


nc 

« * — m- 

m L L 


2(p + q + 1) 

« m m 

-1 + - tan 


-1 


(jt)] = J- <» ♦ p„. q.) , (4) 


where L is the cavity length and p_ and q_ are integers; the notation p , 

m m m 

q is used for labeling frequencies, u * o> . According to the weak 

m mpq 
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coupling approximation, the slowly varying amplitudes A m (t) are complex 
quantities with magnitude A m (t) and phase + m (t)> 

i* m (t) 

V l > = V*> • • <5) 

The superposition (2), involving various frequencies « m , may also contain 
various modal structures P m q m for a given frequency (i.e. degeneracy). 


(a) Particle Equations 

Ve use the guiding center description for the transverse particle 
coordinates 


x * Xg + p cos 0 , 


y = y g + p sin 0 , 


( 6 ) 


p x * p gx ' P 1 sin 9 » p y " p gy + p i cos 9 ’ 

where (x ,y ) and (P^iP^) denote the guiding center position and 
momentum, p^ is the magnitude of the transverse momentum, 0 is the 
gyroangle, p = P_l /0 c * s Larmor radius, « B Q /y is the relativistic 

cyclotron frequency, where fi Q « |e|B o /mc, and the relativistic factor y * 
2 2 1/2 

11 + (pj^/mc) + (p z /mc) j . By averaging the exact Lorentz force 
equations in the vector potential representation over the fast (cyclotron) 
time scale, the slow-time-scale nonlinear relativistic equations of motion 
are cast in the form 


du. 

— ■ E “ n \ Ji<V> sin <*n * V cos Vg • <7a > 

at n 


d0 

dt 


S a> 

o r< n 

“ + E 5" a r 
r n 1 


J l< k n»> 

k n p 


cos (* ♦ 


♦ ) cos k x 
n' n g 


(7b) 
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In Eqs. (7) time has been normalized to w a (S Q /y ) length to 

k Q - - c/w 0 » a n * l e l A n /mc is the normalized radiation amplitude, u is the 

normalized momentum u = p/mc * yv/c and the relativistic factor y = (1 + 
2 2 1/2 

u, + u„ ) . The angle <J» is the relative phase between the m-th mode 

1 z m 

and the particle, 4* = 0 - » t, evolving in time as 
mm 


% 

dt 


= - 


(i) 

E n 
.T a r 


¥ k n p > 

k n p 


cos (4» + 


V cos k n x g 


( 8 ) 


where So> * (w - Q /y) w 1 is the zero-order frequency detuning. The 
nnooo 

prime (') signifies the Bessel function derivative in respect to the 
argument. The evolution of y is found combining Eqs. (7a) to (7b), to 
obtain 


*1 a 
dt 


E - a n J i (k n ,,) sl " <*„ * V cos k n x g • 


(9) 


Two approximations are implicit in Eqs. (7)—(9). The guiding center drift 

has been ignored and the axial relativistic momentum u z is taken as 

constant. They are both justified for a R « 1, kV » 1, since dx^/dt 

2 

dy /dt “ 0(a ) and du /dt ~ 0(a /kV), which are much smaller than du/dt 

& n z n 

de/dt * 0(a n ). 

(b) Field Equations 

Substituting Eq. (2) into Maxwell's equations and ignoring second 
2 “ 2 2 ~ 

order derivatives 3 A /31 << « A , the slow-time evolution of the vector 

n n n 

potentials is given by 
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( 10 ) 


The second term in the right-hand side of Eq. (10) is zero by definition 

for the vacuum eigenmodes defined in Eq. (2b), thus, all the information 

for the slov evolution of the fields is carried in the nonlinear 

transverse current j^. Equation (10) can also be derived from the 

14 

general expression governing coherent wave interaction, 



v 3-1 l _ UU. < 

V g 3x J n c J 


NL 


(ID 


in the limit of standing waves inside the resonator (zero group velocity, 

2 2 

v g = 0) and near vacuum dielectric constant c (weak beam limit: /« o « 

1), (de/d») n a 2u> n . 

Multiplying both sides of Eq. (10) by exp(ia> n t) and taking the fast 

time average over t ~ 2n/(o isolates the n th mode on the left-hand side, 
o o 

Introducing the normalized variables used in Eqs. (7), and letting a n * 
|e|A n /mc = a n exp (i$ n ), ve have 


3a i 4n|e| 

V r > 7T‘ ~ ' <12a) 

3t zu me 

o 

where 

1 _o iw t 

J NL (ri “n> * ~ J 0 dt W 1 -*’ e ' (12b) 

0 


Multiplying both sides by the mode profile U n and integrating over the 
cavity volume lead to 
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( d » . - i 2n|e| f 3 

-" a n- TT- J d r U n (r) V (r? V * (13 > 

'‘St 2Q ' n (0 Z mc 3 V J n NL n 

n on 

3 2 

Above V n = / d r U n is the weighted cavity volume, and the 
phenomenological term « n /Q n vas added to account for resonator losses due 
to boundary effects, where the cavity Q n is defined by 


a 0 iu o 

— a n 2 * - =2- a 2 . (14) 

dt n Q n n 

The current is expressed in terms of the nonlinear distribution 
function f^, 


-2n r “ u. 

j NL ( r »t) = -|e| n b (r) J^d© j^du z j^du^ — cos0 f NL (r,u,t) . (15) 


Letting f NL (r,u,t) = n Q n^r^) f NL^ u l’ U z’ Z;t *’ subs t itu * in g (15) in Eq. 
(13), and expressing r^ m (x,y) and » (u x ,u y ) in the guiding center 
representation Eq. (6), the right-hand side of Eq. (15) becomes 


2n c |e | n 


w 2 me 3 V 
o m 


7^ jSL%L d2 vvv v , g > \ M cos k "’ 


V 6 I 0 dU l I 0 dU i U l ~ COse J l (k n o) WW^'V • 


where again f NL (n« n ) = < 1/T 0 ) J 0 T dt f NL ( r *0 ex P( iw n t >- Expression (16) 
is finally recast in terms of the injected beam current 
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h - i e I c(J zo n o J % I dy E vw ■ 


(17) 


/9 wi . - 

—-2- a = - 2ni < j M . (w ) > . 

L. ) n J NL V n 


3t 2Q 


(18) 


The dimensionless, volume averaged current on the right-hand side of Eq. (18), 


Y 1 n 1 

<j NL < “n> > ‘- 


U C, V T J 
z b n o 0 


o iw t f 
dt e I dr U(r) j NL (r,t) 

J v 


1 Jo 

"z Vn J 


(19a) 


dr U n (r) W'V 


is given, in terms of f NL from Eq. (16), by 


< %<V > ■ 


I. r 


„ dy g n b <x g' y g> h p n (y g ) cos k n X g 


L £ dx gf 

a. V u„ J 0 gJ - 
b n z 

r* r 2 " r" r" > u i 

J J 2 V <2) JJ 6 J n d “a J„ d “i U i J l (k n p) - 0050 £ NL<VV e > Xi “n> ’ 


M n J 0 J 0 * J 0 


3 

where I Q = lell^/mc is the dimensionless beam current (Budker parameter), 

and <j. = JJ dx dy n. (x ,y ) is the effective beam cross section, 
b 8 8 h g g 


(19b) 
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III. NONLINEAR ELECTRON DISTRIBUTION 


The nonlinear electron distribution in phase space f M (u.,u ,0,z,t) 
will be obtained in the form of an asymptotic expansion in powers of the 
small amplitudes a n « 1. Once the form of f^ in the multimode field (2) 
is known, the nonlinear transverse current j ^ is obtained from Eq. (19). 
The evolution of f^ along the particle trajectories Eqs. (7) is given by 
the collisionless kinetic equation 


3f NL 

3f NL 

de 

3f NL 

du z 

3f NL 

du^ 

!5 


4- 

3z 

dt 

39 

+ 1 

dt 

3u z 

+ dt 

3u 


By formally expanding f ^ and ordering in powers of a J , a “ 0(a R ) « 1 


NL 


:(0) 


a f 


( 1 ) 


J f(j) 


one obtains a hierarchy of equations 


:(j) . 


L f VJ = L 1 f 
o 1 


where 


(j-D 


Df 


(j-1) 


I - — 

o " 3t 


u z 3z + 2 c 39 


( 21 ) 


(22a) 


-1 = X 3 " Hi J l^n P) 


i* 3 

i n 

cos k x e - 
n g 


3u, 


1 J,(k p) i\|» 3 

--——— cos k x e n 

2 U± Ir « " « 


k n p 


30J 

(22b) 


Operator L q is the total time derivative L q « d/dt along the unperturbed 
orbits of electrons gyrating in the external magnetic field B q 
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0' - 0 - fi c T , U x ' = U ± , Z 1 = 2 — *~T t * n ' - * n - S« n T ,(23) 

(the characteristic curves of h Q t = 0 in the theory of partial 
differential equations) where 0' a 0(t'), 0 a 0(t) and t = t - t'. 
Replacing the exact trajectories in the right-hand side of (21) by the 
unperturbed orbits (23), (integration along the characteristics), it 
follows that 

f (j) = J dx Df (j_1) [u'(x), 0'(x), z'(t) ] , (24) 

We exploit the fact that the interaction time of an electron inside the 

cavity, of the order of the transit time t^ = W/0 z , is much shorter than 

the characteristic growth time t -1 = d(lna )/dt for the fields. Keeping 

cm * — 

a R constant during the integration (24) yields 


/•v - i(0 - « t.) 

f^^ = E a sin ^ x e ^ 

n n n 8 


(25) 


” dT . h [ z .Vi] e i6 Vj [ J _1*V> !_ + J !< k n p > ±_ L 

} r ’ [ 2i 3u ± (k n p) 2u x 30 J 


:(j-l) 


where T j * ant * we use d 


V'j-l) * e(t j-l ) - Vj-1 ■ 9 j - Vj * *Vj • 

It is now convenient to change variables, introducing the set C Q , Cj, ... 
, C-, defined by 


z u u z u 

C = - , Cj= C- Tj , - , C-= T j = “ ~ ( T i + t 2 + * T i^ 


yW J U Y W 


(26) 
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where dC. = - (u /yV) dx. and C, (x.* 0) ■ C. -. Successive iterations of 
j 2 J J J J- 1 

Eq. (25) yield 


£<J> - l’ d WV jJwV J’.J 0 " dT j L i<V £ 


( 0 ) 


,/rvJ S 


(27) 


'j-l 


J.'WVLwi.J. dt j h<V £ <0) . 

Z 

where is the injected beam distribution, and is written as 

r ^j-l - i(0 - » t.+ v t.) 

J dCj ^(Cj) * E a n « 3 3 cos k_x 


n g 


(28) 


rS- 1 iv n C 

J «, W • 


f J^^p) 3^ J 1 (k n P) 1 3 ) 

2i 3^ (k n p) 2u ± 30 


The quantity v r = (Vy/u z ) S« n , the product of the frequency detuning times 
the transit time through the cavity is the most important parameter for 
the mode coupling coefficients. The integrations, Eqs. (27) and (28), are 
carried out in Appendix A. The expressions for f^ 3 , j*l,2,3 are given by 
(A1)-(A3). For j > 4 the full expressions are getting rather complicated. 
However, one may exploit the existing ordering inside f^ in powers of 
the large parameter h * (kW/y)(u^/u 2 ) » 1. It is then easier to find an 
approximation for f^ to any order in j, Eq. (A6), by keeping the 
dominant contribution. Additional simplifications occur during the 
gyroangle averaging in obtaining the nonlinear current, given by Eqs. (A7) 
and (A8). 
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IV. MODE COUPLING EQUATIONS 


Substituting expression (A7) for the nonlinear current inside the 
right hand side of Eq. (A8) and keeping terms up to fifth order in complex 
amplitudes a ffl = a m exp(i$ m ), the general mode coupling equations assume 
the form 



2ni 

O 0 

A 

v 

dt 

u zo °b ^n 

L 

(i»J»...»m) 

fc< 1 > 

( m;n 

a + C (3) 
a m ^klmjn 

a, a, a + Cjj?-. 
k 1 m ijklmjn 


(29) 


A mode of given frequency « n interacts through all possible frequency 
combinations satisfying the resonant condition 


4 "iJ....in * <“i 1 “j 4 ••• 1 V - “n = 0 ’ < 30 > 

a 

denoted by the notation {i,j,...,m} implies summation over all 

permutations among i,j,...,m (except n) inside Eq. (29), for a given set 

of resonant frequencies (30). A negative index -m will imply -« m , -$ m , 

-v and a * in place of « , ♦ , v and a„ respectively, 
m n m m m n 

The dominant contribution to the coupling coefficients is 

given by 


C 


( 1 ) 

m;n 


2 

V* S 


r* u i 2 r Y w® r 2 df<0) 

■i» J 0 dUl T ) J-l z v " (z) & v " <Cl> ~ 

zo 1 


(31a) 
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,<3) 

■'klm; 


i;n ' (?) 1,4 S kl.;„ “» [‘(Hi.,.' * *ia.,n)] j'du^ (—) 

^ zo 

r“ r 2 ® r^l ^ p ^2 ^ f ^ 

J i z V z > JX v k«l> — J X W — JX V m<S> — ’ 


l l 

(31b) 


5 2 5 

C ijklm;n * ft) ^ S ijklm;n exp [* Xjklmjn 1 + *ijklm;n)] J du ± — 

^ 20 

A 2 v> fi c i vv r-jS yv 5- £■••■ i- £«* w 


df 


( 0 ) 


3u^ -® 3u^ -® 


du, 


where 


V« * J 1 < k m p > h m«> e 


-i\C 

m 


(31c) 


(32) 


♦.. ■(♦.+♦.+...+♦)-♦ is the combined slow phase, and the 

ij...m;n ij mn 

nonlinear filling factor is given by 


L 



$■*(1.$-*... ( i.=y 


2 1 
2 


h i (y g ) h j (y g ) -** h n (y g ) cos k i x g C0S k j X g 


cos k x . 
n g 


(33) 


The results appearing in Eqs. (29)-(33) apply for any spatial 
profile and velocity distribution of the electron beam, as well as any 
radiation modal structure in the open resonator. 
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V. THREE MODE INTERACTION 


Equations (29) are quite general, involving all possible resonant 
combinations of modes. In many cases it has been observed that only a 
small number of modes participates in the final equilibrium with finite 
amplitude. In this section the three mode synchronous interaction, 
satisfying the condition 

\ * V 1 * 2v m • or “k ♦ “l * 2 “m • (34 > 

vill be considered in more detail. This is a special case of the four 
mode interaction, with « m = a> n in Eq. (30). Hence the frequencies w^, o>^ 
are symmetric around « m , \ ® w m ~ du, “l “ “ m + dw » where dw is an y 
frequency interval allowed by the mode spacing in the cavity. The 
possible resonant frequency triads entering are 


W m ' “l - \ > 

- w m - “l + “1 • 

- “m - \ + “k ' 


(35) 


and all 3! permutations of each arrangement, while the possible 
arrangements for are given, setting i * k and j ■ 1, by 


- “m " “1 + “1 * “k + “k 

- w m - “k + “k " “m + w m 


- (i) - w + «_-«_ + « 


“m - “l - \ + “k " “k 

“m - “l - \ + W 1 - “l 

w m ~ “l ' “k + w m ~ w m 


(36) 
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with all the 5! permutations of each. 

There has been a tendency recently to employ the following 

3 

renormalized parameters in gyrotron theory: amplitude F n = a n / Y 0 3 l0 = 
3 2 

E /B 0 , , detuning A = (2/0, )$w , interaction length u = 

n o oi n 10 n 

kW(0 lo 2 /20 zo ) and current I = I Q (8n/V o )(Q/r o 0 lo 3 ) where V Q = Oi/2)k Q 3 LV 2 

3 9 

and I = I b (Ams) je |/m e c x 3 10 . Setting = w^, “3 = “3 ant * 

separating real and imaginary parts, Eqs. (29) are recast in the form 


dF, w. 


TT * ZT F 1 * Rer i F 1 * E , F 1 V ReG l» * F 2 F 3 Re ( G 123 •' J 
at m=i 

i i F, Fj 2 F., 2 R«D llm . i e F 2 2 F 3 Re (D ml23 e- 1 *) , (37.) 

1=1 m=l m=l 


dF„ «, 


^ ^ f 2 - R ' F 2 F 2 * E 3 F 2 V ReG 2„ * F 1 F 2 F 3 R 'l G 231 • J 

i i, F 2 F l 2 F „ 2 ReD 21m * EF m 2 Fj F, F, R.(d. 2j1 .**) , (37b) 
1=1 m=l m=i 


dFj <», 3 


7T F 3 ■ ReF 3 F 3 * 2, F 3 F m 2 ReG lm * F 2 2 F 1 R *( G 312 e ~ **) 
dt zu^ m=i 

* i t F 3 Fj 2 F 2 ReD 3lm . E F „ 2 F 2 2 Fj Re(D m312 e'“) , (37c) 

1=1 m=l m=l 

where T is the complex linear growth rate. The evolution of the slow phase 
+ * +2 + +2 - 2+2 is given by 
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d* 

dt 


SS + Im 


( F 1 


r 3- 


2 r„ 


- F 1 2 Im[G 11 + G 31 - 2G 21 ] - F 3 2 Im[G 13 + G 33 - 2G 23 ) + f^ 1 " 1 ( G 22“ Gj 2 - G 32 ] 
Io ( G 123 e ”) ~ ~~~ Im ( G 312 e “) + F 1 F 3 Im ( G 231 e ' ”) 


F 3 F 2 


F 1 
3 3 


- I, Sj Im Kl« + °31m - 2D 21m) F l 2 


^ “(w”) 


F 1 F 2 


- i# 


Im ( D m312 e if ) * F 1 F 3 Im ( D m231 e 

(37d) 


where 5Q « («^ + « 3 - 2« 2 )/<o o . Note that only one combined slow phase 
appears for each triad. 

The dependence on the gyrotron parameters I, y and A is contained 
in the coupling coefficients. It will be shown now that when y » 
y o 0^ o /2 , it can be eliminated as an independent parameter by proper 
scaling. The computation of r, G, and D, carried out in Appendix B, 
yields 


r n - Vi F <V ’ 
-2 


nml 


J . S 3 5 G < W ’ 

I g S 5 DCv^v.,^) , 


(38) 


n’ V 1' 


.here t . <V> 2(S lo 4/2B z<, 2 > J 1 ' < Vo> and 2 s * 

I(0 lo 2 /8Q)(k o W) 3 (8 io /8 zo ) 3 [J 1 '(k o p o )] 2 . The nonlinear filling factors S r 
S 3 , and Sj are given in (Bll). The quantities T, G, D ,given for a cold 
beam by the integrals (B16)-(B18), depend on the desynchronism parameter 
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\> m alone, the.product of the interaction length y with the detuning A^, 


v = y A 
m m 


(39) 


This suggests the following scaling transformation, valid when « = « and 

n o 

Q n =Q, 


I - 2 I g Q 

P = l F , (40) 

T = t/2Q , 

putting the mode coupling equations in the final reduced form 

^7 ♦ R 1 - i { S l R ' r i F l * S 3 F 1 F . 2 R ' G 1 . * S 3 F 2 2p 3 Re ( G 123 l *) 

* S 5 £j F l F l 2p » 2 “U. - S 5 Be Kl23 e ' 1 ']} • < 4l » 


2 F 2 . i|s lR .r 2 F 2 * S 3 E x F 2 F m 2 Re a 2m * S 3 K.(c 123 e if ) 


dx 


3 3 - . 


* S 5 L m E ., F 2 F l 2 '» 2 R ' D 21» * S 5 J/AVj Re ( D m231 «**)} • < 41b > 


J7 * F 3 ' 2 [ S l ReF 3 F 3 * S 3 V.* ReG 3« * S 3 F 2 2r 1 R '( G 312 e ' **) 


3 3 


3 ' 2 ; i: 


S 5 Sj S 2 VlV Re °31™ * S 5 2 i F « 2 F 2 2 ''i Re (®«.312 •'**)} 


e _1 *U . (41c) 


It follows, from (41), that the normalized start-up current is 1 = 

1/S lV 
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The beam parameters I, 0 ZQ , the interaction length u, the 

cavity Q and the three frequencies combine into only 3 control parameters: 
I, Vj and v ^ (v^ is related to the other tvo frequencies through the 

resonant condition). The interaction length u does not appear explicity, 
but only implicitly through I and v. For simple density profiles, such as 
a pencil or a thin annular beam, the filling factors may also be expressed 
as powers of a single factor s, 

Sj * s^ +1 , s * cos k Q x g h(y g ) , (42) 

with x and y the effective values of x and y . Thus, they could also 
g g g g 

be absorbed in the scaling factor, redefining 

V - s l , 

I' - s 2 I . (43) 

Since I is given analytically in terms of the gyrotron parameters, it 
suffices to compute tables of T, G and D as functions of once, to 
completely know the coupling coefficients for any combination of the 
parameters I, u, 0 ZQ , 0^ o and frequencies w m . Note that the beam spatial 
profile enters only in the nonlinear filling factor S^; in case of 
complicated beam profiles Sj can be computed independently from Eq. (33). 
The discussion so far has been limited to a cold beam. In case of thermal 
effects the general expressions (A7) and (A8) are applicable and the 
control space will increase by two parameters, the pitch angle spread and 
the energy spread. 

Equations (37) and their reduced form, Eqs. (41), are the basic 
result in this paper, and describe three mode interaction up to fifth 
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order in magnitude in the QOG. The diagonal terms G nn describe the 

effects of nonlinear saturation (or self excitation) for each mode. The 

cross-coupling terms G nm describe a phase-independt. . nonlinear damping 

(positive or negative) among different modes. The interaction through the 

G , terms does involve the slow phase and describes mode locking effects, 
nml 

The fifth order terms D nm ^ are necessary to stabilize the system in 
regimes where most of the third order interactions are positive and 
mutually destabilizing. They are also required to account for potential 
amplitude bistability, 1 ^ i.e., the existence of two equilibria of the same 
frequency but of different amplitude. 

Similar results, to third order in amplitude, have been obtained 
for the conventional gyrotron 11 , using a somewhat different approach. The 
present derivation of the coefficients is considerably simpler, since it 
requires a (j+l)-fold integration for the j-th order, compared to the 
(2j+l)-fold integration required in the single particle approach. 11 Mode 
coupling equations, identical in structure, first appeared in the 
treatment of LASER cavities 1 ^. 
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VI. RESULTS AND CONCLUSIONS 


The normalized linear growth rate ReT(v n ) is given in Fig. 2. 
Tables of the third order coupling coefficients G(v n ,v m ) for Gaussian 
profile modes and a cold pencil electron beam have been obtained by 
numerical evaluation of the expressions (B5)-(B7). The contour plots for 
the real and imaginary parts of G nnj , G^ and G 23 J within the regime 0 < 
v n , v m < 3 are shown in Figs. 3a, 3b, and 3c respectively (the 

coefficients G^ are E iven * n terms of G^i G 3 i 2 ^ v l ,v 2 ^ * G 123^ v 3’ v 2^ 

Dotted lines represent regions of negative values for G, corresponding to 

mutually stabilizing influence among interacting modes in that region. 

The mode coupling is destabilizing in the regions with solid lines, 

corresponding to positive G. Note the absence of symmetry in the coupling 

coefficients, where, in general, G„ * G . Contour plots for the fifth 

nm mn 

order coupling coefficients D ( v n * v m * v i> as functions of v m , for 
selected values of v n appear in Figs. 4(a)-4(d). 

According to the final number of participating equilibrium modes, 
Eqs. (41) demonstrate three types of equilibria: (a) single mode, where 
one mode dominates the other two (b) two mode equilibria, where one mode 
is of negligible amplitude and (c) equilibria among three modes of 
comparable amplitudes. The complete set of equilibria among a given set 
of frequencies is given by the zeros of the right-hand side of (42). 
Since F n are defined positive, only solutions with F^, F 2 , F^ > 0 are 

acceptable. It is well known that every stable equilibrium of a nonlinear 
dissipative system, such as the system represented by Eqs. (41), is 
associated with a region of initial conditions in amplitude space (basin 
of attraction) that eventually fall into this equilibrium. Which 
equilibrium the system will choose to settle in cannot be analytically 
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predicted, in the absence of invariants of the motion. This can only be 
done by integration of the set of ordinary differential equations (41). 

The simplest case, a single mode equilibrium for the n-th mode, is 
found from 


d _!n 

dt 


? n ( 


ReT + ReG 


ReD 


0 , 


(44a) 


— = ImT - F 2 ( 
dt n v 


ImG + ImD F 
nn nnn 


n 2 )- 


(44b) 


Equations (44) are obtained by dropping the cross-coupling terms from 
(37); the combined phase # is irrelevant in that case. Acceptable, non¬ 
trivial solutions F n > 0 for the steady-state amplitude should always 
exist, since a single mode saturates through particle trapping. The fifth 

order saturation term D is therefore important on two counts. First, 
nnn 

it is necessary for saturation in regimes where the third order G Rn is 
positive and causes self-excitation instead of suppression. Second, it is 
required to account for amplitude bistability, when two acceptable 

equilibrium values F^ and F^ exist for a given mode, since the third 

1/2 

order nonzero solutions, F r = ± (T /G nn ) , can provide at most one 
acceptable solution F > 0. At steady-state the slow phase $ n varies at a 
constant rate, Eq. (44b); d$ R /dt expresses the nonlinear frequency shift 
from the linear frequency <i> n . 

Equations (37) are used to provide some examples of mode 
competition in a QOG driven by a cold pencil beam in a typical parameter 
regime, with beam current 1^ = 13 A, y q = 1.146, 3^ O /0 ZO = 1* = 5, 
cavity Q « 40,000 and cavity length L = 48cm, corresponding to a frequency 
separation d«/« o * 0.003. In Figs. 4(a)-4(c) we examine mode interaction 
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among the first three frequencies < w 2 < w 3 a ^ ove w 0 > counting only the 
odd resonator modes « n separated by dw/w^. = 0.006. In Fig. 5(a) initially 
? 2 Q » F^ q , F^. The mode F 2 quickly reaches a quasi-steady state, which, 
however, is unstable. The mode F^, which initially grows much slowly than 
its uncoupled linear rate due to nonlinear suppression from F 2 , eventually 
overtakes and suppresses F 2 (mode switching in the cavity). In contrast, 
the steady state reached from F^ q » F 2q , F.j 0 , shown in Fig. 5(b), is 
stable since a large amplitude F^ suppresses F 2 and F^ to negligible 
amplitude. However, the single mode equilibrium of F^ is not accessible 
from the usual start-up initial conditions, vhere all three modes have 
small comparable amplitudes F lo = F 2o “ F 3o « 1. Figure 5(c) shows that 
the final state in this case has F^ as the dominant mode. Note that F 2 
has been completely suppressed in Fig. 5(c), although its linear growth 
rate differs from F 2 by a few percent. Since F^ wins the competition with 
Fj and F 2 , we next examine the coupling of F^ with the next two 
frequencies F^ and F^. It is seen from Fig. 6 that the final steady state 
involves three modes with comparable amplitudes. Note that F^ has been 
excited despite being linearly stable. Figure 7 shows the mode 
competition among the first three modes in a cavity of reduced length, L = 
24cm, and increased frequency separation d w/w q = 0.012 among odd modes. 

In conclusion, the general analytic formalism for multimode 
interaction in the QOG has been developed in this paper. Applications of 
the coupling equations to interactions involving three main modes have 
demonstrated the ability to model nonlinear effects, such as nonlinear 
destabilization, suppression, and mode switching. A new scaling has been 
introduced, during the derivation of the coupling equations, that reduces 
the number of the control parameters. The ratio of operation to start-up 
current I, and the desynchronism parameters for the participating modes 
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become the only free parameters. The inclusion of coupling terms up to 
fifth order in amplitude was shown to stabilize the system in the hard 


excitation regime (where the third order coupling coefficients are 
positive.) The coupling coefficients up to fifth order have been 
tabulated for a cold pencil beam. These results will be generalized in a 
future paper to include thermal spreads in velocity, and arbitrary beam 
current profiles. The set of equations will also be expanded to include 
interactions among a large number of modes, by forming all the possible 
resonant frequency combinations. 
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Appendix A: Computation of the Nonlinear Current 

Substituting operator (28) inside (27) and performing the successive 
iterations yields 



27 



f (3) < v e,c,t> - y E E E >„ \ «J 


l( V V V 1( V V V c 


* ~ * 1 
i a a, < 
n m 1 1 

f 3 

*^1 3 

d M-„" 

J-oo 3u ! J 

d ^2 X. ~ 

3u i 1 


i- £ -n, m, -1 j + £ -n, -m, 1 j + £ -n, -m, -1 J 

* ( V V V i( V V V c 
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.Cl ! 

dC l X - 7" 

3u l J 

d <2 *-m “ 
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dC 3 X 3 


[ -n, m, -1 j + [ -n, -m, 1 J + [ -n, -m, -1 J 


*<V V V i( V V V C 


i 

r c i . 3 


d <2 **“ 
3u i 


d <3*l 


[ -n, -m, -1 ] 


i( V V V i( V V V 1 K 

e e 


d <l X_ n “ 


d <2 X. 


1 


i r^2 

^ L* h Xl 


r i \ 3f<0) 

-n, -m, -1 > - 

1 J J 3u, 


+ cc , 


(A3) 


where the explicit 9, t dependence is contained in = 9 - « n t , 
r , Jc u, x iv ( 

V C ’ U .L>- J l(—J cos k n x g h „<« e ■ X-n <llU l ) * *n <C - u i> 

u o 

zo 

VW - v c -V J iftr) / J i( ! ¥ i ) fir 1 ) - 


In obtaining (A1)-(A3) we used the invariance: + v r C = const., along 

the unperturbed trajectories. The notation ( -n, m, -1 J inside (A2)-(A3) 
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means an integral similar to the preceding fully vritten term except that 
a negative index -n in place of n implies -^ n> -\> n , ^ a R * in place of * n , 

V *-n> V etc ' 

The dominant contribution in each integral comes from the 
derivative 3/3u, acting on v , 


3v « 

- * — h(u,), where 


ntup 


k o Wu l 


Y U„ 


2nV u. 

-- » 

X Y u 2 


1 . 


(A4) 


Therefore, the most important terms inside f^ are those with the 
(3/3u^)^ derivative. Neglecting, for the same reason, terms like 3 y/3u^ - 

if 

3u. /3u, - 1 « 3v /3u, one obtains from (A3) 
x ± m ± 


f < 3 ) ( V e,t,«) - (^) (-^-) E E E Sj a n a.* aj* 


i( V V *1> *<V V V 1 


* 3 

,V 1 3 

dt, v - 

1 “ n 3u, 

J _SB ” U J_ J 

dC 9 v - 

2 " 3u, , 

-® 1 


: ( 0 ) 


[ -n, m, -1 j + [ -n, -m, 1 J + [ -n, -m, -1 J 

iv.C 


+ cc , 


where v (£) * J 
m 


;(¥*) 


h (C) e 

m 


s = cos kx h (y ) . 
m m g m g 


ou l 

(A5) 


One can show by induction, proceeding along the same lines, that 
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where (...) implies all the permutations among the m^' s. 

When the j-th nonlinear piece of is substituted inside (15) 

i8 -i0 

and combined with cos6 » (l/2)(e + e ), the resulting expression 
contains even 9-harmonics, 1, e ± ^ i6 , ...» e ±(j +1 ) i0 f f or j 0 dd f and odd 0- 
harmonics, e ±i0 , e ±3i0 , ...» e ± ^ +1 ^ i ®, for j even. Consequently, coupling 

interactions involving even number of modes vanish completely during the 0- 
integration. The first nonlinear correction is of third order in amplitude, 


the second of fifth, and so on. It should be emphasized that this is a 

consequence of assuming a uniform in 0 initial distribution f^(u); 

however, for a prebunched distribution of the form f^(u,0) = f Q (u) + fj(u) 
16 

e + ... + cc, both even and odd terms survive and the nonlinear effects 


enter to second order in amplitude. Splitting the slowly varying f NL (« n ) 
component off the nonlinear distribution, by multiplying each f^ with 
exp(i« n t) and averaging over the fast time, substituting ^ ns ^ e 
(19), and finally averaging over the gyroangle 0 and over the resonator 
volume, one obtains 
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ff b V n u zo t m i» m 2 ' m 3 * ••• * m 4» n } 
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r" f a r c i a r^2 a rS-i df (0) 

b V2) 1^ d- dC2 ^1.;- «).. d ^’ ~ 


The quantities Aw m v.. , ♦ and S are 

m l ra 2 * * * m j ’n ®j nij m 2* * *n*j;n rn^nij. •*®j »n 

defined in a similar manner as in equations (30), (32) and (33) 

respectively. Of all the possible comoinations [m^n^tm^. . .±ijk J only 
those with (j-l)/2 positive and (j+l)/2 negative signs have survived the 
©-integration and contribute to C^; for example, the term I-l,-m,-n] in 
f^ 3 \ Eq. (A5), has vanished from C^ 3 ^. Also note that the complex 
conjugate terms in (A1)-(A6) drop out during the fast time averaging; they 
contribute to the jj4L^ -w n^ component. 
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Appendix B: Computation of the Three-Mode Coupling Coefficients 

The three-mode coupling coefficients in Eqs. (37) are computed from 
the general expressions in (31), keeping the frequency combinations (35) 
and (36), and yielding, 


f „ ■ K C -”„' 2 

{-m,m,-n} 


.(3) 


h K y r (3) r 

>123 K L -2-23;l’ b 312 
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A cold pencil beam is employed for simplicity, where the velocity 
distribution f^ and the density profile n(Xg,y g ) are S iven b Y 


t (0) 


2 nu 


S( V "io* 5<u z - u zo> > 


lo 


"<YV - 5(x g - f ) S(y g ) 


The calculations are simplified by the following observation involving the 
derivatives of Mu^-u^). If A(u), B(u) are any functions of u satisfying 
d n A/du n « d n B/du n for every h, then 




■ «V 
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and by induction 
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Therefore, moving everything but the exponentials exp(i\> n O in front of the 
velocity integrals, one obtains 
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where 


v n (C) 


iv C 

e MO 


Ve limit the mode coupling among fundamental Gaussian modes where the 
transverse profiles Uqq are 


u 00 (y g ,z) * h 0 ( V h 0 <Z) » h 0 (z) " exp (" ^2 ) » (B9) 

for all frequencies o> n . The normalized volume in that case is given by 


r - r cos k x . 2 2 v re 9 ? 

dx J *>] dz t— jTiTi «!>(- - - k „ k o L “ ■ < bi °> 

y V 

-C 2 iv C 

a. * 1 and v » e e 
d n 

For the pencil beam at x = L/2, and according to Eq. (4), cos k x =0 for 

g n g 

n even and cos k x >1 for n odd. In other words, the electron beam passes 
n g 

through a null of the electric field for modes with an even number of half 
vavelengths betveen the mirrors, or through a maximum, for an odd number of 
half vavelengths. The j-th order nonlinear filling factor Sj, Eq. (33), is 
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then reduced to 
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(Bll) 


Since the frequency separation among the coupled frequencies « n is very 

small, it is important only when the detuning S« n , \> n appears explicitly. 

Thus, one may set k R = k Q , letting Jj'(k n p) * « J Q , = V = 

(n/2)k Q ^LW^ and Q r = Q for all modes. To further simplify (B6)-(B10), one 

may notice that, according to (A4) and for « n = « Q , 

dv kW u. 

— ss q( u .) =- 


du, 


r 


is independent of frequency. Since the S-function causes every derivative 
to be taken at u^ « u 1q , one may employ v r , v m , as independent 
variables and use 
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It is clear from (B16)-(B18) that T, G and D depend only on the initial 
values for the detuning parameters v nQ . Combining (B13)-(B15) with (B16)- 

(B18), converting the amplitudes a R into the normalized F n , a n = F n 

3 3 

(r Q 3^ o ), and absorbing the factors (r 0 f^ 0 ) in the coupling coefficients, 

finally yields 
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and, in general, for the j-th order coupling coefficient, 
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Defining the scaling factor 
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and the new normalized current 
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puts Eqs. (B19)-(B21) in the form given in Eq. (38), 
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Fig. 3 (Continued) — Contour plots for the real parts of third order coupling coefficients (a) G m (b) G, 23 and 
(c) G 23l . Dashed lines signify negative values. 
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(c) 

Fig. 5 (Continued) — Mode competition in the QOG for the first three modes above w„. The time evolution of 
the amplitudes F in (a)-(c) corresponds to different initial conditions. 


50 

























4793/4 DISTRIBUTION LIST 
Updated 2 Aug 89 


Air Force Avionics Laboratory 
AFWAL/AADM-1 

Wright/Patterson AFB, OH 45433 
Attn: Walter Fries 

Air Force Office of Scientific Research 
Bolling AFB 

Washington, D.C. 20332 
Attn: H. Schlossberg 

Air Force Weapons Lab 
Kirkland AFB 

Albuquerque, NM 87117—6008 
Attn: Dr. William Baker 

Dr. Brendon Godfrey 

Bhabha Atomic Research Center 
Laser Division 
Bombay, India 400085 
Attn: T.S. Shirsat 

Columbia University 
520 West 120th Street 

Department of Electrical Engineering 
New York, NY 10027 
Attn: Dr. S.P. Schlesinger 

A. Sen 

Columbia University 
520 West 120th Street 
Department of Applied Physics 
and Nuclear Engineering 
New York, NY 10027 
Attn: T.C. Marshall 

R. Gross 

Cornell University 

School of Applied and Engineering Physics 
Ithica, NY 14853 

Attn: Prof. Hans H. Fleischmann 

John Nation 
r. N. Sudan 

Creol-FEL Research Pavillion 
12424 Research Parkway, Suite 400 
Orlando, FL 32826 
Attn: Dr. Luis R. Elias 

Dr. I. Kimel 


1 copy 


1 copy 


2 copies 
1 copy 


1 copy 


1 copy 
1 copy 


1 copy 
1 copy 


1 copy 
1 copy 
1 copy 


1 copy 
i ..eopy 


53 



Dartmouth College 
18 Wilder, Box 6127 
Hanover, NH 03755 
Attn: Dr. John E. Walsh 

Defense Advanced Research Project Agency/DEO 
1400 Wilson Blvd. 

Arlington, VA 22209 
Attn: Dr. L. Buchanan 

Defense Communications Agency 
Washington, D.C. 20305 
Attn: Dr. Pravin C. Jain 

Assistant for Communications 
Technology 

Defense Nuclear Agency 
Washington, D.C. 20305 
Attn: Mr. J. Farber 

Dr. Leon Wittwer (RAAE) 

Defense Technical Information Center 
Cameron Station 
5010 Duke Street 
Alexandria, VA 22314 

Department of Energy 

Div. of Advanced Energy Projects 

Washington, DC 20545 

Attn: Dr. R. Gajewski 

Department of Energy 
Office of Energy Research 
Washington, D.C. 20545 
Attn: C. Finfgeld/ER-542, GTN 

T.V. George/ER-531, GTN 
D. Crandall/ER-54, GTN 
Dr. David F. Sutter/ER-224, GTN 

Director of Research 
U. S. Naval Academy 
Annapolis, MD 21402-5021 

General Atomics 

13-260 

Box 85608 

San Diego, CA 92138 
ATTN: Dr. J. Doane 

Dr. C. Moeller 


1 copy 


1 copy 


1 copy 


1 copy 
5 copies 


2 copies 


1 copy 


1 copy 
1 copy 
1 copy 
1 copy 


1 copy 


1 copy 
1 copy 


54 







Georgia Tech. EES-EOD 
Baker Building 
Atlanta, GA 30332 

Attn: Dr. James J. Gallagher 1 copy 

Hanscomb Air Force Base 
Stop 21, MA 01731 

Attn: Lt. Rich Nielson/ESD/INK 1 copy 

Hughes Aircraft Co. 

Electron Dynamics Division 
3100 West Lomita Boulevard 
Torrance, CA 90509 
Attn: J. Christiansen 

J. Tancredi 

Hughes Research Laboratory 
3011 Malibu Canyon Road 
Malibu, CA 90265 
Attn: Dr. R. Harvey 

Dr. R.W. Schumacher 

KMS Fusion, Inc. 

3941 Research Park Dr. 

P.O. Box 1567 
Ann Arbor, MI 48106 
Attn: S.B. Segall 

Lawrence Berkeley Laboratory 
University of California 
1 Cyclotron road 
Berkeley, CA 94720 
Attn: Dr. A.M. Sessler 

Lawrence Livermore National Laboratory 


P.O. Box 808 
Livermore, CA 94550 

Attn: Dr. D. Prosnitz 1 copy 

Dr. T.J. Orzechowski 1 copy 

Dr. J. Chase 1 copy 

Dr. W.A. Barletta 1 copy 

Dr. D. L. Birx 1 copy 

Dr. E.T. Scharlemann 1 copy 

Litton Electron Devices 
960 Industrial Road 
San Carlos, CA 94070 

Attn: Library 1 copy 


1 copy 
1 copy 


1 copy 
1 copy 


1 copy 


1 copy 


55 





Los Alamos National Scientific Laboratroy 
P.0. Box 1663, MSJ 564 
Los Alamos, NM 87545 

Attn: Dr. Brian Newman 1 

Los Alamos Scientific Laboratory 
P.0. Box 1663, AT5-827 
Los Alamos, NM 87545 

Attn: Dr. T.J.T. Kwan 1 

Dr. L. Thode 1 

Dr. C. Brau 1 

Dr. R. R. Bartsch 1 

Massachusetts Institute of Technology 
Department of Physics 
Cambridge, MA 02139 

Attn: Dr. G. Bekefi/36-213 1 

Dr. M. Porkolab/NW 36-213 1 

Dr. R. Davidson/NW 16-206 1 

Dr. A. Bers/NW 38-260 1 

Dr. K. Kreischer 1 

Dr. B. Danby 1 

Dr. G.L. Johnston 1 

Massachusetts Institute of Technology 
167 Albany St., N.W. 16-200 
Cambridge, MA 02139 

Attn: Dr. R. Temkin/NW 14-4107 1 

Mission Research Corporation 
8560 Cinderbed Road, Suite 700 
Newington, VA 22122 

Attn: Dr. M. Bollen 1 

Dr. J. Pasour 1 

Naval Research Laboratory 
Addressee: Attn: Name/Code 

Code 0124 - ONR 1 

Code 1000 - Commanding Officer 1 

Code 1001 - T. Coffey 1 

Code 1003.9A - Computer Resources Architect 1 

Code 1005 - Head, Office of Mgt. & Admin. 1 

Code 1005.1 - Deputy Head, Off. of Mgt. & Admin. 1 

Code 1005.6 - Head,_Directives Staff - 1 

Code 1200 - Capt. R. W. Michaux 1 

Code 1201 - Deputy Head, Command Support Division 1 

Code 1220 - Security 1 

Code 2000 - J. Brown 1 

Code 2604 - NRL Historian 1 

Code 2628 - TID Distribution 22 

Code 2634 - Cindy Sims 1 

Code 3000 - R. Doak 1 


copy 


copy 

copy 

copy 

copy 


copy 

copy 

copy 

copy 

copy 

copy 

copy 


copy 


copy 

copy 


copy 

copy 

copy 

copy 

copy 

copy 

copy 

copy 

copy 

copy 

copy 

copy 

copies 

copy 

copy 


56 





Code 4000 - W. Ellis 
Code 4600 - D. Nagel 
Code 4700 - S. Ossakow 
Code 4780.1 - W. Ali 
Code 4790 - Branch Office 
Code 4794/4793 - Library 
Code 4793 - W. Black 
Code 4770 - G. Cooperstein 
Code 4794 - A. Fliflet 
Code 4793 - S. Gold 
Code 4790 - C. Hui 
Code 4795 - C. Kapetanakos 
Code 4793 - A. Kinkead 
Code 4790 - Y. Lau 
Code 4707 - W. Manheimer 
Code 4794 - M. Rhinewine 
Code 4790 - P. Sprangle 
Code 5700 - L. Cosby 
Code 6840 - S. Ahn 
Code 6840 - A. Ganguly 
Code 6840 - R. Parker 
Code 6840 - N. Vanderplaats 
Code 6850 - L. Whicker 
Code 6875 - R. Wagner 

Naval Sea Systems Command 
Department of the Navy 
Washington, D.C. 20362 
Attn: Commander, PMS 405-300 

Northrop Corporation 
Defense Systems Division 
600 Hicks Rd. 

Rolling Meadows, IL 60008 
Attn: Dr. Gunter Dohler 

Oak Ridge National Laboratory 

P.O. Box Y 

Mail Stop 3 

Building 9201-2 

Oak Ridge, TN 37830 

Attn: Dr. A. England 

Office of Naval Research 
800 N. Quincy Street 
Arlington, VA 22217 
Attn: Dr. C. Roberson 

Office of Naval Research 

1012 W 36th Street, Childs Way Bldg. 

Los Angeles, CA 90089-1022 

Attn: Dr. R. Behringer 


1 copy 
1 copy 
26 copies 
1 copy 
4 copies 
20 copies 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 


1 copy 


1 copy 


1 copy 


1 copy 


1 Copy 


57 






Optical Sciences Center 
University of Arizona 
Tucson, AZ 85721 

Attn: Dr. Willis E. Lamb, Jr. 1 copy 

Physical Sciences, Inc. 

635 Slaters Lane #G101 
Alexandria, VA 22314-1112 

ATTN: Dr. M.E. Read 1 copy 

Physics International 
2700 Merced Street 
San Leandro, CA 94577 

Attn: Dr. J. Benford 1 copy 

Princeton Plasma 
Plasma Physics Laboratory 
James Forrestal Campus 
P.O. Box 451 


Princeton, NJ 08544 

Attn: Dr. H. Hsuan 

2 

copies 

Dr. D. Ignat 

1 

copy 

Dr. H. Furth 

1 

copy 

Dr. P. Efthimion 

1 

copy 

Dr. F. Perkins 

1 

copy 

Raytheon Company 

Microwave Power Tube Division 

Foundry Avenue 

Waltham, MA 02154 

Attn: N. Dionne 

1 

copy 

Sandia National Laboratories 

ORG. 1231, P.O. Box 5800 

Albuquerque, NM 87185 

Attn: Dr. Thomas P. Wright 

1 

copy 

Mr. J.E. Powell 

1 

copy 

Dr. J. Hoffman 

1 

copy 

Dr. W.P. Ballard 

1 

copy 

Dr. C. Clark 

1 

copy 

Science Applications, Inc. 

1710 Goodridge Dr. 

McLean, VA 22102 

Attn: Adam Drobot 

1 

copy 

P. Vitello 

1 

copy 

D. Bacon 

1 

copy 

C. Menyuk 

1 

copy 

Science Research Laboratory 

15 Ward Street 

Somerville, MA 02143 

Attn: Dr. R. Shefer 

1 

copy 


58 






SPAWAR 

Washington, D.C. 20363 
Attn: E. Warden, Code PDE 106-3113 

Capt. Fontana, PMW 145 

Spectra Technologies 
2755 Northup Way 
Bellevue, WA 98004 
Attn: Dr. J.M. Slater 

Stanford University 

Dept, of Electrical Engineering 

Stanford, CA 94305 

Attn: Dr. J. Feinstein 

Stanford University 
High Energy Physics Laboratory 
Stanford, CA 94305 
Attn: Dr. T.I. Smith 

TRW, Inc. 

One Space Park 
Redondo Beach, CA 90278 
Attn: Dr. H. Boehmer 

Dr. T. Romisser 
Dr. z. Guiragossian 

University of California 
Physics Department 
Irvine, CA 92717 
Attn: Dr. G. Benford 

Dr. N. Rostoker 

University of California 
Department of Physics 
Los Angeles, 90024 

Attn: Dr. A.T. Lin 

Dr. N. Luhmann 
Dr. D. McDermott 

University of Maryland 
Department of Electrical Engineering 
College Park, MD 20742 
Attn: Dr. V. L. Granatstein 

Dr. w. W. Destler 


1 copy 
1 copy 


1 copy 


1 copy 


1 copy 


1 copy 
1 copy 
1 copy 


1 copy 
1 copy 


1 copy 
1 copy 
1 copy 


1 copy 
1 copy 


59 




University of Maryland 
Laboratory for Plasma and Fusion 
Energy Studies 
College Park, MD 20742 
Attn: Dr. Tom Antonsen 

Dr. John Finn 
Dr. Jhan Varyan Heilman 
Dr. W. Lawson 
Dr. Baruch Levush 
Dr. Edward Ott 
Dr. M. Reiser 

University of New Mexico 
Department of Physics and Astronomy 
800 Yale Blvd, N.E. 

Albuquerque, NM 87131 
Attn: Dr. Gerald T. Moore 

University of Tennessee 
Dept, of Electrical Engr. 

Knoxville, TN 37916 
Attn: Dr. I. Alexeff 

University of Utah 

Department of Electrical Engineering 
3053 Merrill Engineering Bldg. 

Salt Lakfe City, UT 84112 
Attnt \|Dr. Larry Barnett 
'/Dr. J. Mark Baird 

U. S. Army 
Harry Diamond Labs 
2800 Powder Mill Road 
Adelphi, MD 20783-1145 
Attn: Dr. Howard Brandt 

Dr. Edward Brown 

Dr. Stuart Graybill 

Dr. A. Kehs 
Dr. J. Silverstein 

Varian Associates 
611 Hansen Way 
Palo Alto, CA 94303 
Attn: Dr. H. Huey 

Dr. H. Jory 
Dr. Kevin Felch 
Dr. R. Pendleton 
Dr. A. Salop 

Varian Eimac San Carlos Division 
301 Industrial Way 
San Carlos, CA 94070 
Attn: C. Marshall Loring 


1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 
1 copy 


1 copy 


1 copy 


1 copy 
1 copy 


1 copy 
1 copy 
1 copy 
1 copy 
1 copy 


1 copy 
1 copy 
1 copy 
1 copy 
1 copy 


1 copy 





Yale University 
Applied Physics 
Madison Laboratory 
P.0. Box 2159 
Yale..Station 
New Haven, CN 06520 
Attn: Dr, I. Bernstein 

Ken Busby 


Director of Research 
U.S. Naval Academy 
Annapolis, MD 21402 
(2 copies? 


Do NOT make labels for 
Records-(01 cy) 


Naval Research Laboratory 
Washington, DC 20375-5000 
Code 2630 
Timothy Calderuiood 






